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ABSTRACT 

a Context. Circumstellar disc evolution is paramount for the understanding of planet formation. The gas in protoplanetary discs large 

program (GASPS) aims at determining the circumstellar gas and solid mass around ~ 250 pre-main-sequence Herbig Ae and T Tauri 
_^ stars. 

Oj' Aims. We aim to understand the origin and nature of the circumstellar matter orbiting 5 1 Oph, a young (< 1 Myr) luminous B9.5 star. 

I ' Methods. We obtained continuum and line observations with the PACS instrument on board the Herschel Space Observatory and 

O ■ continuum data at 1.2 mm with the IRAM 30m telescope. The spectral energy distribution and line fluxes were modelled using the 

physico-chemo radiative transfer code ProDiMo to constrain the gas and solid mass of the disc around 5 1 Oph. The disc vertical 

C/J ' hydrostatic structure was computed self-consistently together with the gas thermal balance. 

C^ ' Results. We detected a strong emission by atomic oxygen [OI] at 63 microns using the Herschel Space Observatory. The [OI] 

emission at 145 microns, the [CII] emission at 158 microns, the high-y CO emissions, and the warm water emissions were not 
detected. Continuum emission was detected at 1.2 mm. The continuum from the near- to the far-infrared and the [01] emission are 
well explained by the emission from a compact (iJout=10-15 AU) hydrostatic disc model with a gas mass of 5x10"* Mq, 100 times 
that of the solid mass. However, this model fails to match the continuum millimeter flux, which hints at a cold outer disc with a mass 
in solids of ~ 10"* Mq or free-free emission from a photoevaporative disc wind. This outer disc can either be devoid of gas and/or is to 
cold to emit in the [O I] line. A very flat extended disc model (i?ou,=400 AU) with a fixed vertical structure and dust settling matches 
all photometric points and most of the [O I] flux. 

Conclusions. The observations can be explained by an extended flat disc where dust grains have settled. However, a flat gas disc 
^^ , cannot be reproduced by hydrostatic disc models. The low mass of the 51 Oph inner disc in gas and dust may be explained either by 

f^ . the fast dissipation of an initial massive disc or by a very small initial disc mass. 
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1 . Introduction the small dust grains result from the collisions of planetesimals. 

lZ , The dust discs in the latter category are geometrically flat and 

C^ . The large variety of extrasolar planetary systems requires a ^j^^ ^^^ ^j^^^^ ^^^y. ^^^^^^ -^ ^j^^ ^^^..jr ^j^^ ^^^^ 

■ ■ theoretical and observational understanding of protoplanetary ^^^ ^^ ^^^ ^^^^^ categories seems rapid because there are very 

discs where the exoplanets are formed. Observing protoplane- ^^^ ^^ ^^^ ^^.^^jj^^ transition discs. The definition of a transi- 

tary discs at different stages of evolution makes it possible to ^-^^ ^-^^ -^ ^^^ commonly accepted yet, but discs with a lack of 

test the various models of disc dissipation, planet-formation, and ^^^._^^ ^^^^^^ ^^^^^^^ ^^ ^^ -^^^^ j^^j^ ^^^^^^ ^ ^ j^^^^^ ^^^ 

planet-migration. Discs are characterized by their gas and solid ^j^^^ ^^^^ ^ f^^.j^ ^^^^^^ ^^^ ^^ considered as transitional 

contents and their geometry. The two main categories of discs ^-^^^ ^j^^ -^^ ^^^^^^ ^^^^^^^^ B9.5IV-B9.5V star 5 1 Oph 

are the young massive and flanng ones (Md.c=10--10-^ M^ ^^ j^^ j^j^ ^59^^ H^ 158(543^ ^^q ^^^^jq^ Ir^S 17283-2355) 

m gas and dust, typical size of 100-500 AU), charactenzed by j^^^^^^ at -131 pc is a rare example of a pre-main-sequence star 

a strong IR excess in their spectral energy distribution (SED); ^j^^^ -^ surrounded by a disc where most or all of the primor- 

and the older debns discs, where there is no gas left and where ^^^j ^^^ ^^^ p^^^^^j^ dissipated. The SED of the 51 Oph disc 

~~r~r, TT^ TTTT , .... suggests a new category of disc. It does show a near-IR excess, 

Herschel is an ESA space observatory with science instruments , , . Vr- 1 or-i^ • 1 • 1 

■ , , , c 1 J n ■■ 1 T *■ t -*• J fu ■ but the main property of the SED is a drop in the continuum 

provided by European-led Principal Investigator consortia and with im- f f j f 

portant participation from NASA ^^^ beyond 20 microns. No observations existed beyond 60 mi- 
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Fig. 1. Herschel spectra towards 51 Oph. The [OI] line is spectrally unresolved. The line peak wavelengths are shown. 



crons. The absence of strong far-IR continuum dust emission 
points to a small amount of dust grains and a non-flaring disc. 
The 10-microns feature is consistent with emission by amor- 
phous silicates and, there are almost no features that can be 
attributed to polyc y clic aromatic hydr ocarbon (PAH) emissions 
(iKeller et al.ll2008h . [Stark et al.l (l2009l) proposed a two-opticafly 
thin disc model to explain the combination of MIDI-VLT, Keck 
nulling, and Spitzer-IRS data. Their outer dust disc extends to 
1200 AU and is flat. 

The evolutionary status of 5 1 Oph is debated. The low dust 
mass fav ours the idea that 5 1 Oph is a young debris disc as re- 
ported in lStark et al.l (l2009l) . A debris disc would have a lower 
gas-to-dust mass ratio than a young disc. One of the key ques- 
tions is therefore how much gas (atomic + molecular) is still 
present in the 51 Oph disc? The total atomic + molecular gas 
mass of the 5 1 Oph disc has never been determined. In addition 
to the mass, the nature of the gas in the 5 1 Oph disc is also dis- 
puted. Is the gas atomic or molecular? A predominately atomic 
disc would favour the scenario where most of the gas comes from 
the evaporation of infalling comets without primordial gas left, 
whereas significant amount of molecular gas (H2) as traced by 
rotational emissions by CO and its isotopologues would indicate 
that the gas is a leftover from planet-formation and/or disc dissi- 
pation. Absorptionstudieshi the ultraviolet with the FUSE tele- 
scope dGradv & Silvisll 19931: iLecavelier Pes Etangs et"anil997t 



iRoberge et al.ll2002h suggest infalling atomic gas and no molec- 
ular gas, similar to the case of the yS Pic debris disc. On the 
other hand, spectra in the infrared show ro-vibrational emis- 
sion lines from molecules such as CO, CO2, and H2O (e.g. 
Ivan den Ancker et al.l200lh . The presence of molecules suggests 
that the gas may be dense and warm enough so that efficient gas- 
phase neutral-neutral reactions can produce molecules such as 
water. Neutral carbon (C I) has been detected in absorption with 
the FUSE telescope, and the C I/dust column density ratio is low, 
similar to the value found in the line-of-sight of /? Pic. Hot car- 
bon monoxide (CO bandhead emission) has been observed by 
several groups and was assigned to the emission from a hot gas 
in a disc seen close to edge-on (iThi et al.ll2005l: iBerthoud et al.l 
2007 ; Tatulli et al. 2008) , but no rotational emission from cold 
(r=10-30 K) CO in the sub-millimeter range is available. In this 
paper, we present new observations in the millimeter and far-IR 
obtained with the IRAM 30m telescope and the Herschel-PACS 
Space Observatory in Sec. |2| & E] The 51 Oph Herschel data 
is part of the HerbigAe and debris disc survey from the GASPS 
large open time programme (P.I. Dent; Dent et al., in press). The 
new observations (continuum and line) combined with archival 
data are modelled to derive an estimate of the disc dust and gas 
masses in Sec. |4| In the discussion, we consider the nature and 
the different possible origins of the 51 Oph disc (Sec.|6]l. 
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2. Observations and data reduction 

The source 5 1 Oph was observed for contin uum emission at 1.2 
mm using the MAMB02 bolometer array (iKrevsaet al.lll998h 
on the IRAM 30m telescope at Pico Veleta, Spain. Observations 
were conducted during the bolometer pool in 2008 (Nov. 13-18, 
23). Zenith opacity for our observations was typically ~0.2-0.3. 
Observations were first carried out to a target 1 cr sensitivity of 
1 mJy (20 minutes on-source) and then repeated on two subse- 
quent nights to reach 5 sigma detection. A total of 60 minutes 
on source was achieved, in an ON-OFF pattern of 1 minute on 
target followed by 1 minute on sky, with a throw of 32". Flux 
calibration was carried out using Mars, Saturn, or Jupiter (de- 
pending on availability), and local pointing and secondary flux 
calibration was carried out using IRAS 16293-2422B. The data 
were reduced using the facility reduction software, MOPSICQ. 

We obtained photometric and line observations with the 
PACS instrument dPoghtsch et al. 20101) on board the Herschel 
Space Observatom (JPilbratt et al.l l2010h . We obtained PACS 
scan map observations at 70, 100, and 160 microns. At 100 
microns, 51 Oph was observed at two angles (70 and 110 de- 
grees) to improve the noise suppression, while at 70 microns, 
we only observed one scan direction. Both the 70 and 100 mi- 
crons scans include the red band at 160 microns. The PACS pho- 
tometric data were reduced with the mini scanmap pipeline in 
HIPE version 8.1.0, (calTree version 32), and we constructed 
mosaics of the two scan maps at 100 and 160 microns. The con- 
tinuum flux was extracted using an aperture of 12" (sky between 
30 and 35") at 70 microns, of 15" (sky between 33 and 38") 
at 100 microns, of 20" (sky between 40 and 50") at 160 mi- 
crons. We obtained fluxes of 873+23 mJy at 70, 347+10 mJy 
at 100 and 90±6 mJy at 160 microns (see Table [TJ. The cali- 
bration errors haven been included. The spectroscopic data were 
reduced with the official version 8.0.1 of th e Herschel Interactive 
Processing Environment (HIPE: IOtg|20 10b . using standard tasks 
provided in HIPE. These include bad pixel flagging; chop on/off 
subtraction; spectral response function division; rebinning with 
oversample = 2 and upsample =1, corresponding to the native 
resolution of the instrument; spectral flatfielding and finally av- 
eraging of the two nod positions. To conserve the best signal and 
to not introduce additional noise, we only extracted the central 
spaxel and corrected for the flux loss with an aperture correction 
provided by the PACS instrument team. We extracted the fluxes 
of the detected lines using a Gaussian fit to the emission lines 
with a first-order polynomial to the continuum, using the RMS 
on the continuum (excluding the line) to derive a 1 cr error on 
the line by integrating a Gaussian with height equal to the con- 
tinuum RMS and width of the instrumental FWHM. In case of a 
non-detection, we give a 3 cr upper limit, also calculated on the 
continuum RMS. The measured line fluxes are listed in Table [1] 

Pipeline-processed FEROSpESO 2.2m high-resolution (R~ 
45000) optical spectra (4000-9000 A) were obtained from the 
Advanced Data Products ESO archiveQ. Details of the FEROS 
data reduction pipeline can be found in the ESO HARPS and 
FEROS pipeline web-site Q The spectra are displayed at the ob- 
served wavelength and no correction was performed for the ra- 
dial velocity or motion of the earth around the sun. Finally, we 
complemented the new data with archival photometry. 



51 Oph Ha 



' http://www.iram.es/IRAMES/mainWiki/CookbookMopsic 

^ The Fiber-fed Extended Range Optical Spectrograph, Kaufer, A. et 
al. 1999, The Messenger 95, 8 

^ http://archive.eso.org/ 

■* http://www.eso.org/sci/facilities/lasilla/instruments/harps/tools/ 
software.html/ 
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Fig. 2. Normalized optical spectrum of 5 1 Oph around the Ha 
emission obtained by FEROS. 
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Fig. 3. Modified blackbody fits to the 51 Oph SED. The error 
bars on the photometric points are 3 cr. The orange line is the 
Spitzer-IRS spectrum taken from the NASA archive. The best fit 
withy6=-0.3 is unphysical. 



3. Observational results 

We detected the spatially unresolved (on a scale of 10") emis- 
sion at 1.2mm with a flux of 5 ± 1 mJy (1 cr error) towards 51 
Oph using the MAMBO bolometer at the IRAM 30m telescope, 
consistent with the up per limit of 42 mJy obtained at 1 .1mm by 
ISvlvester et al.l (Il996l) . The continuum was detected by PACS in 
all settings. The total SED from the optical to the millimeter is 
shown in Fig.fTl lFaiardo-Acostaet al.](ll993h detected the 10 //m 
silicate featur e at low spectral resolution. The feature was later 
confirmed bv Ivan den Ancker et al.l (|2001|) . The millimeter flux 
translates into a dust grain (with radii < 1mm) mass of (1-2) x 
10"^ Mq assuming rdust= 50 K and the standard dust opacity 
A-y ^ 2.0(v/vo)^g''cm^ wh ere vq = 230.769 GHz andyS = 0.6 
dBeckwith & Sargentll99lh . This derived value is an upper limit 
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Fig. 4. Disc gas density, gas-to-dust ratio, SED, and localization of the continuum emission for the fixed vertical structure disc 
model. The disc is gas-rich above the mid-plane and gas-poor in the mid-plane. For clarity a couple of observed photometric points 
are not shown in the SED. The lower-right panel shows the cumulative continuum emission for 15% and 85% for each wavelength. 
The continuum emission for wavelengths greater than 500/zm comes from the outer disc {R >100 AU). 



to the dust mass since we cannot estimate the contribution from 
free-free emission at 1 .2mm due to the lack of data at longer 
wavelengths. Alternatively, pebbles with size greater than 1 cm 
will emit as blackbodies in the millimeter range. Assuming a 
gas-to-dust mass ratio of 100, the disc gas mass is (1-2) x 10""* 
Mo. 

We complemented the SED with archival photometric data 
and the Spitzer-IRS spectrum (see Fig. [3]l. We fitted the SED 
with a Kurucz stellar model in addition to a modified blackbody 
of the form 



VFy = By{T){A/Ao) 



-P-i 



(1) 



where Bv{T) is the Planck function at the temperature T . In 
the formula, we assumed that the grains will emit with k - kq 



g"' cm" for A < Aq and with Ky - Ko{A/A()y^ g"' cm^ for A > Aq. 
The value of Aq is set by the criterion Aq > In < a >, where 
< fl > is the mean grain radius. We adopted /lo = 100 fim (< a>- 
16 yum). The fits to the SED for three values of yS (/?=!, 0, -0.3) 
are shown in Fig. [3] A v alue for fi of 1 is typical for grains found 
in protoplanetary discs (lDrainell2006l) . A pure blackbody emis- 
sion has y6=0. Alternatively, a grain size powerlaw distribution 
with an index P = 3 w ill result in a quasi-constant k with wave- 
length (lDrainell2006l) . The model with the standard value for 
protoplanetary discs fi—\ fails to match the millimeter flux. A 
pure blackbody emission also underpredicts the millimeter flux, 
but the prediction is still consistent with the observations within 
3 cr. The best fit is obtained with y8=-0.3. A value for/? of -0.3 
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Fig. 5. Disc gas density structure from hydrostatic models 1 and 2 (7?out=2 AU, Mdisc=5xlO"* M©) between 2 and 15 AU. The 
left panel corresponds to the model with dust settling (model 1), the right panel to the model without dust settling (model 2). The 
dust temperature contours are overplotted in black. The red contours correspond to the Stokes numbers (5, =0.1, 1, and 10 contour 
values). 
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Fig. 6. Disc dust density structure from hydrostatic models 1 and 2 (/?out=2 AU, Mdisc=5xlO"^ Mq) between 2 and 15 AU. The left 
panel corresponds to the model with dust settling (model 1), the right panel to the model without dust settling (model 2). The dust 
temperature contours are overplotted in black. The Stokes parameters are overplotted in red (S ,-Q.\, 1, and 10 contour values). 



is unphysical and points to the emission by an extra component, 
either a cold outer disc/envelope or free-free emission. 

In addition, far-IR spectral line observations were obtained 
using the Herschel Space Observatory. The only detected line is 
the [OI] emission at 63 microns. The [OI] emission is confined 
to the central spaxels. This line is always the str ongest in the 
two [ O I] and [C II] surveys of the GASPS project dMeeus et al.l 
I2OI2I) . No high-y CO or water lines were detected. The line flux 
and three sigma upper limits are summarized in Table [1] Finally, 



we show the archival Ha line profile in Fig.|2] The profile shows 
a double-peaked structure on the top of the stellar Ha absorption. 
The origin of this profile is discussed in Sec.|4] 

4. SED and gas modelling 

We modelled the continuum and gas emission with the codes 
MCFOST and ProDiMo. MCFOST is a continuum and NLTE 
line Monte Carlo radiative transfer code that can handle 2D 
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Fig. 7. Grain mass absorption (/Cab.s in red dotted- line), scattering 
(/fsca in black dashed-line), and extinction (/c in black full line) 
opacity for the adopted dust grain properties listed in Table |2] 
The slope in the upper-right corner shows a hypothetical opacity 
with/f ~ /t '. 
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Fig. 8. Fit to the SED for model 1 with and without dust settling. 
The four far-IR photometric points are Herschel-PACS observa- 
tions. The 1 .2 mm photometry was obtained with MAMBO at the 
IRAMSOm telescope. 



and 3D geometries (iPinte et al.ll2006i l2009l) . ProDiMo is a ra- 
diative chemo-physical code that solves the continuum radia- 
tive transfer, the gas heating and cooling balance, the gas chem- 
istry, and the disc vertical hydrostatic stmcture self-consistently 
with the gas pressure (Woitke et al. 2009). New features o f the 
codes are described in lKamp et alj ( l2010l) . lThi et al.l (1201 Ih . and 
IWoitke et al.l (l201lh . The 2D continuum radiative transfer mod- 




r[AU] 

Fig. 9. Gas temperature structure for model 1 . The gas tempera- 
ture is higher than 250 K in the inner compact disc. 




r LAUJ 

Fig. 10. [OI] emission from the disc for model 1. The upper 
panel shows the continuum at 63 microns and the line verti- 
cal optical depth. The middle panel shows the cumulative [OI] 
63 micron emission from the inner to the outer disc. The lower 
panel exhibits the atomic oxygen density distribution in the disc. 
The contour encloses the 70% x 70% emission area. 



ule in ProDiMo an d the code MCFO ST has been benchmarked 
against other codes dPinte et al.ll2009h . 

The modelling strategy is twofold. First we assumed a fixed 
vertical structure disc with dust settling. The disc e xtends from 
1.5 AU to 400 AU as suggested by the model of IStarketaP 
(l2009l) . We employed an automatic fit procedure based on the 
Simplex optimization algorithm. The free parameters of the 
model (model a) are the total disc gas-ndust mass, the surface 
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Table 1. Lines and photometric data obtained by Herschel- 
PACS. The errors are Icr and the upper limits are 3cr. The cali- 
bration error adds an extra ~ 10% uncertainty. 



T able 2. Common dis c parameters. The parameters are defined 
in lWoitkeet"aD (120091) . 
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Cent, flux 


Line flux 
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(Jy) 


(10-"* Wm-2) 


[0I]3Pi ^^^Pj 


63.183 


1.09 + 0.19 


53.3 ±2.5 


[OI]3Po^3p^ 


145.525 


0.11 ±0.10 


< 5.1 


[CII]2P3/2^^Pl/2 


157.74 


0.07 ±0.14 


< 6.8 


CO y=36 ^ 35 


72.85 


0.71 ± 0.20 


<11.6 


CO y=33 ^ 32 


76.36 


0.52 ± 0.28 


<16.1 


CO y=29 ^ 28 


90.16 


0.50 ±0.19 


< 8.0 


COy=18^ 17 


144.78 


0.12 ±0.07 


< 4.2 


OH 1/2 - 3/2 


79.11 


0.54 ±0.19 


<11.4 


OH 1/2 - 3/2 


79.18 


0.50 ±0.19 


<11.3 


O-H2O 818 - 7o7 


63.32 


1.08 ±0.15 
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O-H2O 7o7 - 616 
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0.71 ±0.23 
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P-H.0 322 -2ii 
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P-H2O 4n - 32^ 


144.51 
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P-H2O331 -4o4 


158.31 


0.09 ±0.15 


< 7.3 
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< 5.7 


CH+ y = 5 - 4 


72.14 


0.64 ± 0.3 
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Herschel-PACS continuum data 
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0.347 ± 0.010 






160.0 


0.090 ± 0.006 





density profile index, the gas-to-dust ratio, the disc gas scale- 
height, and the flaring index. In model b, the free parameters are 
the same except for the gas-to-dust ratio, which is fixed at 100. 
The dust scale height including settling is 



H'{a,r) - //(r)max[l,fl/fls] 



-j/2 



(2) 



where H{r) is the gas scale height, and s and a^ are two fixed pa- 
rameters: s is the power-law index and a^ is the minimum radius 
for a grain t o be aff'ected by settl ing. The settling law is described 
in details in lWoitke et al.l (|2009|) and has previously been used in 
modelling the HD 163296 disc ( iTilUng et alj|2012l) . We adopted 
a minimum dust radius for settling as of 0.5 //m and s-\. The 
disc density and gas-to-dust ratio structure are shown in Fig.|4] 



Second, we relaxed the assumption of fixed disc vertical struc- 
ture by letting the ProDiMo code calculate the gas scale-height 
that would be consistent with the local gas pressure as deter- 
mined by the gas heating and cooling balance. 

The disc dust and gas parameters are gathered in Table |2] 
The stellar parameters are taken from iMontesinos et al.l (l2009h . 
The input stellar spectrum is taken from the Kurucz database 
of models. The disc is seen close to edge-on between 80 °-90° 
(iThi et al.ll200l iBerthoud et al.ll2007t iTatulli et al.ll2008h . The 
optical photometry shows that the star is seen without strong 
circumstellar extinction, therefore the dust disc has to be flat. 
On the other hand, the gas disc should be more extended verti- 
cally s uch that absorption lines can be de t ected (iGradv & Silvisl 
Il993t iLecaveUer Pes Etangs et al.l IT997t iRoberge et al.l l2002h . 
Therefore we tested disc models where the gas and dust are well 
mixed and where the dust grains have settled towards the mid- 
plane. 



Stellar mass 


M. 


3.8 Mo 


Stellar luminosity 


L, 


260 Lo 


Effective temperature 


Teff 


10000 K 


Distance 


d 


130 pc 


Disc inclination 


i 


80 ° (0 °face-on) 


Disc inner radius 


Rm 


1.5 or 2.0 AU 


Disc outer radius 


^Ollt 


15 or 400 AU 


Vert, column density index 


e 


1 


Inner rim soft edge 




on 


Gas-to-dust mass ratio 


6 


100 


Dust grain material mass density 


Pdust 


3.5 g cm^^ 


Minimum dust particle size 


^^nin 


0.1 ;/m 


Maximum dust particle size 


^max 


5000 iim 


Composition 




ISM silicate 


Dust size distribution power law 


P 


3.5 


H2 cosmic ray ionization rate 


(CR 


1.7x10 -"s-' 


ISM UV field w.r.t. Draine field 


X 


1.0 


Abundance of PAHs relative to ISM 


/pah 


10-3 


a viscosity parameter 


a 


0.0 


Turbulence width 


8v 


0.15kms-i 



Total disc mass 


Md,,sc 


5x10-^-1x10-* Mo 


Dust settling power law 


s 


0, 1.0 


Minimum grain radius for settling 


fls 


0.5 yum 



The accretion rate of the 5 1 Oph disc is di s puted . 
Ass uming a magn e tosphe ric origin, iGarcia Lopez et all (|2006|) 
and iBrittain et al.l (l2007h estimated a gas accretion rate of 
(l-2)xl0~^ M(7) y r-' fr om the hydrogen Bry line. However, 
iMendigutia et alj ( 1201 ll) found that the Wa 10% peak width 
for 51 Oph is much lower than the expected minimum value 
required by the magnetospheric accretion model (Fig. |2). 
Therefore, they argued that the inner disc of 5 1 Oph is decou- 
pled from the star and that the mass accretio n rate of 10'^ Mp 
yT_Ms an upper limit. In a subsequent paper, IMendigutia et al.l 



( I2OI2I) derived an empirical relation between the disc mass /Wdis 
in Mq and the disc mass accretion /Wacc in Mq yr"' valid for both 
T Tauri and Herbig Ae stars 



log/W,ec = l.l(+0.3)xlog/V/disc -5.0(+0.5). 



(3) 



Using this relation and upper limit to the gas disc mass of 2x 1 0""* 
Mq derived from the 1 .2mm flux assuming a gas-to-dust mass 
ratio of 100, we obtained an upper limit to the mass accretion 
rate of ~ 8.5x10-'" M© yr-'. In this paper, we modelled the disc 
around 51 Oph as passive {a viscosity parameter is in Table|2). 
The Ho- emission most likely does not trace magnetospheric ac- 
cretion. The origin of the double-pe aked Hg ernission profile 
can be ascribed to a stellar/disc wind (iKraus et al.ll2008t) . A low 
turbulent width of 0.15 km s-' was assumed. 

The inner disc radius is consistent with the location of the 
sublimation radius for a grain size distribution that follows a 
power-law with index 3.5 (a,nin=0.1 /im and flniax=2000 jum). 
The dust grains are composed of pure amorphous silicate 
(Astronomical silicates, Laor& Draine 1993). The choice of 
amorphous silicate is motivated by the absence of strong crys- 
talline silicate features in the infrared (iKeller et al.ll2008l) . The 
dust opacities computed using the Mie theory for compact spher- 
ical particles are shown in Fig. |7] Except from 10 and 100 jum, 
the opacity is dominated by the scattering term. 

Even a perfect fit to the SED alone does not allow one to 
constrain the disc size, which is a key parameter in many phys- 
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ical phenomena : dust coagulation rate, dust-gas interaction, etc 
(lArmitagell2010l) . One of the main uncertainties is the disc outer 
radius. We adopted two values for the outer disc: a small disc 
with /?out=15 AU and an large disc with /?out=400 AU. The high 
va lue for the outer ra dius at 400 AU was motivated by the model 
of IStark et alj (|2009|) while a compact disc is suggested by the 
drop in the continuum flux at wavelengths beyond 20 micron. 
We also modelled various discs with a sub-standard gas-to-dust 
ratio (models 9 to 16). 

The OI emission at 63.18 yum probes only relatively warm 
and dense gases because of its upper energy level at 227 K and 
a critical density of ~4.7xl0^ cm"^. In the absence of detected 
lines that probe cold gases, we cannot constrain the gaseous con- 
tent beyond 10-20 AU well. 



Table 3. Best-fit model parameters for the fixed vertical disc 
structure model. 



Total disc gas mass 


Meas+dust (Mo) 


1.1x10-* 


Gas-to-dust mass ratio 




29, 100 


Abundance of PAHs relative to ISM 


/pah 


10-3 


Surface density profile index 


e 


0.4 


Reference gas scale height 


HoiAU) 


0.14 


Reference radius 


Rr.f (AU) 


5.0 


Flaring index 


y 


0.55 



One of the major unknown in the gas line modelling is the 
amount of PAHs in the disc. The Spitzer-IRS spectr um in Fig. [3] 
does not show the typ ical IR features of PAHs dKeller et alJ 
l2008HStark eraDl2009l) . We ran the models with the PAH abun- 
dance depleted with respect to the interstellar value of 3 x 10"^ 
by a factor lO"-'. 

The second source of uncertainties are the elemental abun- 
dances relative to H-nuclei. In young protoplanetary discs where 
the gas comes for the collapse of a prestellar core, the elemen- 
tal abundances are the same as for molecular clouds. For debris 
discs, the gas comes from the evaporation of icy planetesimals 
and the gas is enriched in elements heavier than hydrogen and 
helium. The evolutionary state of 5 1 Oph is disputed. The pres- 
ence of significant amount of hot CO as evidenced by the band- 
head emission, of water vapour, and of carbon dioxide support 
the hypothesis that the oxygen and carbon elemental abundance 
with respect to hydrogen are close to the molecular cloud values. 

The chemical network was limited to 7 1 gas and ice species 
including formation of H2 on grain surfaces and fast neutral re- 
action involving vibrationally excited H2. Most of the reaction 
rates are taken from the UMIST database (WoodaU et al. 2007) . 



We adopted a standard value of 1 .7 x 1 
ionization rate. 



s for the cosmic -ray 



The code takes cooling by water lines from the near-IR to 
the far-IR into account. T he coUisional data w ere taken from the 
Leiden-Lambda database dSchoier et al.ll2005h . The original ref- 
erences for the experimental or theoretical rates are given in ap- 
pendix [A] The observed [O I] flux is best reproduced by a com- 
pact (15 AU) disc with masse of SxlO"* Mq. The disc gas-to- 
dust mass ratio is consistent with the interstellar value of 100. 
All other modelled fluxes are below the 3 sigma upper limits. 

We ran two groups of vertical hydrostatic models: one with 
a compact disc and one with a large disc. The parameters of the 
models are listed in Table |4] 



5. Model results 

We first present the results from the Simplex automatic fitting us- 
ing the fixed vertical structure disc model. The fitting procedure 
leads to an extremely flat disc for the gas {Hlr-2%). The dust 
disc is even flatter because of settling (upper panels of Fig. |4]i. 
The best model has a gas mass of IxlO""* M© and a dust mass of 
1x10"^ M0, similar to the values derived from the simple anal- 
ysis of the millimeter emission. The best model fits to the SED 
and the continuum emission cumulative plots of 15% and 85% 
are shown in the lower panel of Fig. |4] The cumulative plots 
show that the continuum emission up to ~100 yum comes from 
within the first 20 AU of the disc. Most of the [O I] emission is 
confined within 20 AU as well. 

For the vertical hydrostatic discs, none of the models, either 
compact (7?out=15 AU) or extended (/?out=400 AU), managed to 
simultaneously fit the SED from near-IR to the millimeter flux 
density. The best simultaneous fit to the SED from the near- to 
the far-IR, excluding the data point at 1.2 mm and [OI] flux, is 
obtained for a compact disc (model 1 and 2) with a disc gas mass 
of 5x10-'' Mq and a dust mass of 5x10-^ Mq. The fits to the 
SED with and without dust settling do not fit the observed con- 
tinuum emission at 1 .2mm. The models underpredict the 1.2 mm 
flux because of the small dust mass in the models. Models 7 and 
8, both with dust mass of 10-^ M© as suggested by the simple 
analysis of the 1.2 mm flux, overpredict all observed fluxes at 
wavelengths shorter than 1.2 mm. The gas scale-heights in the 
hydrostatic models are much higher than the value suggested by 
the best model with a fixed vertical structure disc. As a result, 
even with dust settling, the hydrostatic models of 400 AU discs 
produce far too strong far-IR continuum fluxes. 

In the absence of spatially resolved image at 1 .2 mm, we can 
constrain neither the location of the cold dust grains nor their 
geometry p recisely. However, an extended flat disc has been ad- 
vocated by IStark et al.l (l2009l) . who fitted the SED and mid-IR 
interferometry data simultaneously. An alternative explanation 
for the 1.2 mm emission is stellar free-free emission, which can 
also explain the Ha and Bry emission. 

Hydrostatic disc models with dust settling provide a better 
fit of the mid- and far-IR photometric points while the model 
without dust settling performs better in the near-IR but less well 
in the far-IR. This conclusion on the settling applies only for 
the assumption that the grains are composed of pure compact 
spherical amorphous silicate grains. Other grain compositions 
and shapes are possible. Grains with amorphous carbon inclu- 
sions emit more efficiently in the near-infrared. The [OI] flux 
from the unsettle-model is closer to the observed value but both 
flux estimates are well within the observational uncertainties. We 
tested the effect of well-mixed gas and dust (no settling) model 
and all other parameters the same. The fit to the SED of the well- 
mixed model is slightly worse (oveipredicting the flux in the far- 
IR) than in the model with settling and the [OI] emissions are 
comparable, albeit different (Table|4]i. 

The disc vertical hydrostatic structure shown in Fig. |5] is 
complex but clearly drops beyond 8 AU due a drop in the gas 
temperature in both the settled and well-mixed models. The dust 
density structure follows that of the gas structure (Fig. |6). The 
dust grains are warmer farther out in the settled-dust model, 
which is expected since the grains are on average smaller in the 
disc upper layers. We also show the gas temperature in Fig. |9] 
which is much higher than the dust temperature in the disc at- 
mosphere. 

The [OI] flux from models 1 and 2 slightly unpredicts the ob- 
served flux. The [OI] flux emission area is shown in Fig. [10] The 
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Table 4. Model parameters and line fluxes with 3 cr upper limits. 



Model 


"out 


Mdis,. 


gas-to-dust 


settling 


[OI]63yum 


[Ol] 145 yum 


[Cll] 157 //m 




(AU) 


(Mo) 


ratio 




(10-"* Wm 


-') 


















Observations 














53.3 + 


2.5 


<5.1 


<6.8 














Fixed vertical str-ucture model 


a 


400 


1.0x10-* 


29 


yes 


35.8 




1.3 


0.17 


b 


400 


LOxlO-* 


100 


yes 


65.4 




3.6 


0.02 














Vertical hydrostatic models 


1 


15 


5x10-" 


100 


yes 


47.2 




2.8 


0.02 


2 


15 


5x10-" 


100 


no 


49.4 




2.9 


0.03 


3 


400 


5x10-" 


100 


yes 


22.2 




0.6 


1.2 


4 


400 


5xlO-« 


100 


no 


22.2 




0.6 


1.2 


5 


400 


5x10-' 


100 


yes 


65.3 




2.8 


1.1 


6 


400 


5x10-5 


100 


no 


64.9 




2.7 


1.0 


7 


400 


IxlO-"* 


100 


yes 


94.6 




4.2 


1.2 


8 


400 


IxlO-"* 


100 


no 


93.7 




4.2 


1.2 


9 


15 


5x10-' 


10 


yes 


38.2 




1.1 


0.03 


10 


15 


5x10-' 


10 


no 


40.8 




1.2 


0.03 


11 


15 


5x10-** 


1 


yes 


5.3 




0.1 


0.005 


12 


15 


5x10-'* 


1 


no 


5.3 




0.1 


0.005 


13 


400 


5x10-5 


10 


yes 


144.7 




4.2 


9.8 


14 


400 


5x10-5 


10 


no 


134.6 




3.9 


8.4 


15 


400 


5x10-" 


1 


yes 


39.5 




1.0 


2.5 


16 


400 


5x10-" 


1 


no 


39.5 




1.0 


2.5 



[O I] emission at 63 micron arises from the disc midplane within 
the first 10-15 AU. The [OI] line does probably not probe the 
same circumstellar region as the 1 .2 mm continuum emission. 
It is clear that the [O I] line cannot be used to probe cold gas 
emission. 

In all models, atomic oxygen is by far the main oxygen car- 
rier in the disc, followed by CO. The disc midplane is molec- 
ular and CO reaches its maximum abundance of 1.3x10-"*. In 
our models the [O I] emission traces the disc molecular gas. The 
disc around 51 Oph is mostly molecular (H2, CO, etc) and poor 
in C^ apart from the upper disc atmospheres. The resulting C II 
fluxes are low, consistent with the observed upper limit. There is 
an intermediate layer where neutral atomic carbon is abundant. 
The oxygen 63 micron line is emitted within 15 AU of the star 
regardless of the physical size of the disc. 

The main gas coolant in the [O I] emitting region are the CO 
ro-vibrational lines. The [OI] line probes a tiny fraction of the 
total disc gas mass only. In the absence of PAHs the main heating 
is the CO IR-pumping and background heating by atoms and 
ions (Fell , Sill, CII, CI) in the first 5 AU and dust photoelectric 
on dust grains beyond 5 AU. Water is a very minor component, 
and its ro-vibrational lines do not contribute to the gas cooling. 

A larger disc (7?out=400 AU) with the same gas mass as the 
best compact model predicts the continuum and [O I] line flux. 
The modelling confirms that we cannot constrain the actual disc 
outer radius. However, the location of the OI emission is con- 
fined to the first 10-15 AU for both the compact and the more 
extended disc. The more massive large discs (models 7 and 8) do 
not match the SED and overpredict the [O I] flux. A few models 
with sub-standard gas-to-dust mass ratio have [O I] flux close to 
the observed values. However, like all other hydrostatic model, 
the fits to the SED are poor. 



6. Discussion 

What is the origin of the gas and dust grains in the inner disc? 
Do the grains come from the coagulation of interstellar grains, 
or are they the results of collisions between planetesimals scat- 
tered from the outer belt reservoir by unseen companions? Is 
the gas th e remnant of an initial m assive disc? 51 Oph is young: 
~0.7 Myr jMontesinos et alj2009h . while the youngest stars with 
established debris discs are about 10 Myr old (j6 Pic has an age 
of 12 Myr). Could the disc around 51 Oph be considered as a 
debris disc at its age? 

The best fixed vertical structure model is a flat gas disc 
whereas the gas scale height derived from the hydrostatic mod- 
els is much more extended vertically. The comparison between 
the fixed vertical structure disc model and the hydrostatic mod- 
els suggests that there may be a mechanism to flatten the gas 
disc around 5 1 Oph. The vertical hydrostatic structure is derived 
from balancing the thermal pressure and the gravity forces. 

The disc average gas-to-dust mass ratio is ~ 30-100 for the 
best model, consistent with the standard interstell ar value and 
with the ratio in other H erbig Ae discs: HD 169142 (iMeeus et al.l 
l2010h and HD 163296 dTilling et alJl2012h . A gas-to-dust mass 
ratio of 100 is much high er than the ratio in the debris discs 
(e.g., iLebreton et al.l l2012h . where the dust grains come from 
shattering of planetesimals and the gas is absent. Nevertheless, 
the 51 Oph disc gas mass (5x10-^ M©) is much lower than the 
gas mass of the aforementioned Herbig Ae discs, which can reach 
0.01 Mq. If the disc around 51 Oph is primordial, as suggested 
by its young age, then either the initial disc had a low mass or the 
initial massive disc has been actively dissipating either by pho- 
toevaporation, accretion onto the star, and/or planet formation. 
As discussed above, the current mass accretion rate onto 5 1 Oph 
is difficult to assess, not to mention the accretion rate history. In 
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discs where the gas is replenished by planetesimal evaporations, 
the gas-to-dust ratio would be closer to unity. Hydrodynamic 
simulations show that dust grains can gather into eccentric rings 
as observed in near-IR ima ges of debris discs wit hout the need 
of the presence of a planet dLvra & Kuchnerll2012l) . 

The high luminosity of the star {L^- 260-300 Lq) compared 
with other HerbigAe stars (L«= 20-40 L©) may explain a fast 
disc photoevaporation or increases in the efficiency of radiative 
pressure on grains resulting in a radially extended disc. Disc 
photo-evaporation occurs when hot gases have enough energy to 
escape the star-i-disc gravitational potential entraining the dust 
grains that are dynamically coupled to the gas. The upper disc 
atmosphere outside the critical ra dius of ~6.8 AU for 51 Oph 
is unbound and flows as a wind dHollenbach et alj|1994l) . The 
mass-loss rate through disc evaporative wind is Mw - 4 x 10"'" 
Mq yr"' assuming a stellar ionizing flux of <1) = 10"*^ ioniz- 
ing photons per second, valid for a low-mass solar-lu minosity 
star (lAlexander & Armitagel2009llPascucci et al.l201 ll) . Scaling 
to the luminosity of 51 Oph, the wind mass-loss rate reaches 
1.2x10"^ yr '. A detailed disc-wind model is needed to test 
whether the mass-loss rate of ~10"^ yr"' can explain the ob- 
served Ha and Bry line fluxes and profiles. In 0.7 Myr, the disc 
around 51 Oph would have lost 0.1 Mq of gas. If photoevapo- 
ration is occurring at a constant rate, the initial disc gas mass 
around 51 Oph would have been ~0.1 Mq compared with the 
current stellar mass of 2.8 Mq. Since the inner disc is warm, the 
initial disc would have been gravitationally stable. The current 
disc within 10-15 AU will quickly disappear at this high pho- 
toevaporation rate. A mass-loss rate of ~10"^ yr"' is very high 
compared with the mass in in the disc (~ 5x10"^ M©) and may 
be overestimated. A detailed disc-wind model is warranted to 
better constrain the photoevaporative mass-loss rate. 

The dust mass in the small compact disc around 5 1 Oph de- 
rived from the simultaneous fit to the SED (near- to far-IR pho- 
tometry data) using the hydrostatic disc model is small at 5x10"** 
Mq (0.017 Me), two orders of magnitude lower than if all the 
1 .2 mm flux were attributed to dust emissi on. However, this dust 
mass i s still ten times the value derived bv lvan den Ancker et al.l 
( 1200 lb Their value was derived from a fit to photometric points 
at wavelengths shorter than 100 microns only. Our derived dust 
mass is five times higher than the mass found by iTatuUi et al.l 
(I2008h . 

The discrepancies can also be attributed to differences in dust 
opacities. The continuum emission at 63 micron is optically thin, 
thus it probes the amount of cool dust in the disc. 

The best models with a fixed vertical structure or with a disc 
in vertical hydrostatic equilibrium suggest that the gas-to-dust 
mass ratio is ~ 30-100, but are the gas and dust well-mixed or 
have the dust grains settled to the midpl ane? We used the mode l 
results to compute the Stokes parameter (IWeidenschillingll 977h . 

_ Q.K agmd 

In CsPg 

where Q/f - {GM^IR^y^- is the Keplerian angular velocity at 
radius R, ad the grain radius in cm, nid is the volume mass den- 
sity of a dust grain, which is ~ 1 g cm~^ for a silicate grain, c, the 
sound speed in cm s"' , andpj = «hX a.m.u. the gas mass density 
in g cm"-' (a.m.u. is the atomic mass unit in grams). The Stokes 
parameter varies throughout the disc. It measures the level of 
dust grain coupling with the gas via drag force. For 5/ »1, the 
gas drag is inefficient and the dust grains quickly settle towards 
the disc midplane. For 5, <*:!, the gas and dust grains are well 
mixed and are collocated. The contours for logjQ the Stokes pa- 
rameters (values of -10, 1, and 10) in the disc models with and 



without settling for the mean grain size are plotted over the den- 
sity distribution in Fig. |5] and over the dust density distribution 
in Fig. |6] The Stokes parameter is small throughout the disc in 
our models because the gas is dense («h > 10*' cm"-') and warm 
(T^,, > 250 K). 

Fits to the SED and gas lines do not strongly constrain the 
presence of dust settling. However, the SED together with VLT- 
MIDI and Keck n ulling data support models with dust settling 
( IStark et al.ll2009l) . A more detailed modelling is needed to com- 
pute the exact dust settling, coagulation , and migration, which 
is bey ond the scope of this paper (e.g. iBarriere-Fouchet et al.l 
l2005h . 

Interestingly, dust settling does not much affect the gas ver- 
tical hydrostatic structure for two reasons First, photoelectric- 
heating by dust grains is not the only heating process in the disc 
layers located above the midplane. Heating by PAHs, which are 
small enough to be coupled to the gas, is important, but their 
abundance is very low. Second, sub-micron grains, which are 
more strongly coupled to the gas and remain in the disc atmo- 
sphere, provide most of the photoelectrons. Likewise, the gas 
line collisional excitations do not change much in our model of 
the 5 1 Oph disc whether the dust grains have settled or not. 

If we can attribute the millimeter emission to a cold outer 
disc instead of free-free emission and the near- to mid-IR con- 
tinuum emission and the [OI] to a warm compact inner disc, 
the circumstellar environment of 5 1 Oph resembles a young and 
massive counterpart to main-sequence stars with a warm exozo- 
diacal dust ring and a cold outer planetesimal belt (iDefrere et al.l 
I2OI2I) . The main difference between 51 Oph and fi Pic is that 
the dust grains in an exozodiacal ring around fi Pic have to be 
replenished, while part the dust grains around 5 1 Oph may still 
be primordial. 

7. Conclusions 

We detected [OI] far-IR and 1.2 mm continuum emissions to- 
wards 5 1 Oph. The [O I] probes the warm gas in a compact in- 
ner circumstellar disc while the 1 .2 mm emission suggests a cold 
outer disc. An extended flat disc with dust settling can explain 
the continuum and line data consistent with previous model. 
Hydrostatic disc models predict discs that are too extended in 
the vertical direction for the gas and the dust. The failure of 
the hydrostatic models to explain the observations suggests that 
mechanism(s) other than thermal pressure and gravity forces are 
at work to keep the gas and dust disc flat. The gas-to-dust mass 
ratio is consistent with the standard value of 100. Future observa- 
tions in the (sub)-millimeter with ALMA will help constraining 
the cold gas and dust around 5 1 Oph. 

Acknowledgements. We thank the Agence Nationale pour la Recherche (con- 
tracts ANR-07-BLAN-0221, ANR-2010-JCJC-0504-01, and ANR-2010-JCJC- 
0501-01), CNES, PNPS of CNRS/INSU, France, and support from the 
Millenium Science Initiative (Chilean Ministry of Economy), through grant 
"Nucleus P10-022-F" for support. We acknowledge funding from the EU 
FP7-2011 under Grant Agreement nr 284405 (PERG06-GA-2009-256513). 
Computations presented in this paper were perfonned at the Service Commun 
de Calcul Intensif de I'Observatoire de Grenoble (SCCI) on the super-computer 
Fostino funded by ANR. C. Eiroa, G. Meeus, and B. Montesinos are partly 
supported by Spanish grant AYA 2011-26202. We finally thank the referee for 
his/her comments that helped improving the manuscript. 

References 

Abrahamsson, E., Krems, R. V., & Dalgamo, A. 2007, ApJ, 654, 1 171 
Alexander, R. D. & Armitage, P J. 2009, ApJ, 704, 989 

Armitage, P. J. 2010, Astrophysics of Planet Formation (Cambridge University 
Press) 



10 



W.-F. Thi: The disc around 51 Oph 



Barber, R. J., Tennyson, J., Harris, G. J., & Tolchenov, R. N. 2006, MNRAS, 

368, 1087 
Barriere-Fouchet, L., Gonzalez, J.-R, Murray, J. R., Humble, R. J., & Maddison, 

S. T 2005, A&A, 443, 185 
Beckwith, S. V. W. & Sargent, A. 1. 1991, ApJ, 381, 250 
Bell, K. L., Berrington, K. A., & Thomas, M. R. J. 1998, MNRAS, 293, L83 
Berthoud, M. G.. Keller, L. D., Herter, T L.. Richter, M. J., & Whelan, D. G. 

2007, ApJ, 660, 461 
Brittain, S. D., Simon, T, Najita, J. R., & Rettig, T W. 2007, ApJ, 659, 685 
Chambaud, G., Levy, B., Millie, R, et al. 1980, Journal of Physics B Atomic 

Molecular Physics, 13, 4205 
Daniel, E, Dubemet, M.-L., & Grosjean, A. 2011, A&A, 536, A76 
Defrere, D., Lebreton, J., Le Bouquin, J.-B., et al. 2012, A&A, 546, L9 
Draine, B. T 2006, ApJ, 636, 1114 
Dubemet, M.-L. & Grosjean, A. 2002, A&A, 390, 793 
Fajardo-Acosta, S. B., Telesco, C. M., & Knacke, R. F. 1993, ApJ, 417, L33 
Faure, A., Grimier, N., Ceccarelli, C., et al. 2007, A&A, 472, 1029 
Faure, A., Gorfinkiel, J. D., & Tennyson, J. 2004, MNRAS, 347, 323 
Flower, D. R. 2001, Journal of Physics B Atomic Molecular Physics, 34, 2731 
Flower, D. R. & Launay. J. M. 1977, Journal of Physics B Atomic Molecular 

Physics, 10, 3673 
Garcia Lopez, R., Natta, A., Testi, L., & Habart, E. 2006, A&A, 459, 837 
Grady, C. A. & Silvis, J. M. S. 1993, ApJ, 402, L61 

Hammami, K., Owono Owono, L. C., & Stauber, P 2009, A&A, 507, 1083 
Hollenbach, D., Johnstone, D., Lizano, S., & Shu, E 1994, ApJ, 428, 654 
Jankowski, P & Szalewicz, K. 2005, J. Chem. Phys., 123, 104301 
Jaquet, R., Staemmler, V., Smith, M. D., & Flower, D. R. 1992, Journal of 

Physics B Atomic Molecular Physics, 25, 285 
Kamp, 1., Tilling, I., Woitke, P, Thi, W.-E, & Hogerheijde, M. 2010, A&A, 510, 

A18 
Keller, L. D., Sloan, G. C., Eon-est, W. J., et al. 2008, ApJ, 684, 41 1 
Kraus, S., Hofmann, K.-H., Benisty, M., et al. 2008, A&A, 489, 1157 
Kreysa, E., Gemuend, H.-P, Gromke, J., et al. 1998, in Society of Photo-Optical 

Instrumentation Engineers (SPIE) Conference Series, Vol. 3357, Society 

of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, ed. 

T.G.Phillips, 319-325 
Laor, A. & Draine, B. T 1993, ApJ, 402, 441 
Launay, J.-M. & Roueff, E. 1977, Journal of Physics B Atomic Molecular 

Physics, 10, 879 
Lebreton, J., Augereau, J.-C, Thi, W.-E, et al. 2012, A&A, 539, A17 
Lecavelier Des Etangs, A., Vidal-Madjar, A., Backman, D. E.. et al. 1997, A&A, 

321, L39 
Lim, A. J., Rabadan, I., & Tennyson, J. 1999, MNRAS, 306, 473 
Lyra, W. & Kuchner, M. J. 2012, ArXiv e-prints 

Meeus, G., Montesinos, B., Mendigutia, 1., et al. 2012, A&A, 544, A78 
Meeus, G., Pinte, C, Woitke, R, et al. 2010, A&A, 518, L124 
Mendigutia, 1., Calvet, N.. Montesinos, B., et al. 2011, A&A, 535, A99 
Mendigutia, 1., Mora, A., Montesinos, B., et al. 2012, A&A, 543, A59 
Montesinos, B., Eiroa, C, Mora, A., & Merin, B. 2009, A&A, 495, 901 
MuUer, H. S. P 2010, A&A, 514, L6+ 
Offer, A. R., van Hemert, M. C, & van Dishoeck, E. E 1994, J. Chem. Phys., 

100, 362 
Ott, S. 2010, in Astronomical Society of the Pacific Conference Series, Vol. 434, 

Astronomical Data Analysis Software and Systems XIX, ed. Y. Mizumoto, 

K.-l. Morita, & M. Ohishi, 139 
Pascucci, 1., Sterzik, M., Alexander, R. D., et al. 201 1, ApJ, 736, 13 
Pilbratt, G. L., Riedinger, J. R., Passvogel. T, et al. 2010, A&A, 518, LI 
Pinte, C, Harries, T J., Min, M., et al. 2009, A&A, 498, 967 
Pinte, C, Menard, E, Duchene, G., & Bastien, P 2006, A&A, 459, 797 
Poglitsch, A., Waelkens, C, Geis, N., et al. 2010, A&A, 518, L2 
Roberge, A., Feldman, P. D., Lecavelier des Etangs, A., et al. 2002, ApJ, 568, 

343 
Schoier, E L., van der Tak, E E S., van Dishoeck, E. E, & Black, J. H. 2005, 

A&A, 432, 369 
Stark, C. C, Kuchner, M. J., Traub, W. A., et al. 2009, ApJ, 703, 1 188 
Sylvester, R. J., Skinner, C. J., Barlow, M. J., & Mannings, V. 1996, MNRAS, 

279, 915 
Tatulli, E., Malbet, E, Menai'd, E, et al. 2008, A&A, 489, 1151 
Thi, W.-E, van Dalen, B., Bik, A., & Waters, L. B. E M. 2005, A&A, 430, L61 
Thi, W.-E, Woitke, P, & Kamp, I. 201 1, MNRAS, 412, 711 
Tilling, I., Woitke, R, Meeus, G., et al. 2012, A&A, 538, A20 
Turpin, E, Stoecklin, T, & Voronin, A. 2010, A&A, 511, A28+ 
van den Ancker, M. E., Meeus, G., Cami, J., Waters, L. B. E M., & Waelkens, 

C. 2001, A&A, 369, L17 
Weidenschilling, S. J. 1977, MNRAS. 180, 57 
Wemli, M., Valiron, P, Faure, A., et al. 2006, A&A, 446, 367 
Wilson, N. J. & Bell, K. L. 2002, MNRAS, 337, 1027 
Woitke, P, Kamp, I., & Thi, W.-E 2009, A&A, 501, 383 



Woitke, P, Riaz, B., Duchene, G., et al. 201 1, A&A, 534, A44 

Woodall, J., Agundez, M., Markwick-Kemper, A. J., & Millar, T J. 2007, A&A, 

466, 1197 
Yang, B., Stancil, R C, Balakrishnan, N., & Forrey, R. C. 2006, J. Chem. Phys., 

124, 104304 



Appendix A: Collisional data 

The original articles for the line frequen cies, E i nstein 
coefScients, and colli s ional rates are CO (Flower 20011 
Jankowski & Szalewiczl l2005t lYang et al.1 12006; .Wernli et al. 
20061) . H 7O (iBarberetaU 120061: iDubernet & Grosiean 
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